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SPECTROPOLARIMETRY AND MODELING OF THE ECLIPSING T TAURI STAR KH 15D 
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ABSTRACT 

KH 15D is a strongly variable T Tauri star in the young star cluster NGC 2264 that shows a decrease in 
flux of 3.5 magnitudes lasting for 18 days and repeating every 48 days. The eclipsing material is likely 
due to orbiting dust or rocky bodies in a partial ring or warped disk that periodically occults the star. 
We measured the polarized spectrum in and out of eclipse at the Keck and Palomar observatories. 
Outside of the eclipse, the star exhibited low polarization consistent with zero. During eclipse, the 
polarization increased dramatically to ~ 2% across the optical spectrum, while the spectrum had the 
same continuum shape as outside of eclipse and exhibited emission lines of much larger equivalent 
width, as previously seen. From the data, we conclude that (a) the scattering region is uneclipsed; (b) 
the scattering is nearly achromatic; (c) the star is likely completely eclipsed so that the flux during 
eclipse is entirely due to scattered light, a conclusion also argued for by the shape of the ingress and 
egress. We argue that the scattering is not due to electrons, but may be due to large dust grains 
of size ~ 10/.tm, similar to the interplanetary grains which scatter the zodiacal light. We construct 
a warped-disk model with an extended dusty atmosphere which reproduces the main features of the 
lightcurve, namely (a) a gradual decrease before ingress due to extinction in the atmosphere (similar 
for egress); (b) a sharper decrease within ingress due to the optically-thick base of the atmosphere; 
(c) a polarized flux during eclipse which is 0.1% of the total flux outside of eclipse, which requires 
no fine-tuning of the model. The inclination of the warp is set by the duration of ingress and egress 
assuming that the warp is located at the Keplerian radius, and the inclination of the observer is then 
determined by the duration of the eclipse. 

Subject headings: eclipses — polarization — stars: pre-main sequence — stars: variables: other — 
stars: individual (KH 15D) 



1. INTRODUCTION 

T Tauri stars are thought to be young stars that 
are st ill accreting gas through a disk. iKearns fc Herbstl 
(1998) discovered T Tauri star KH 15D (K7 V, d = 760 
pc), which shows an eclipse for 1/3 of its period, im- 
plying that it is obscured by faint circumstellar material 
rather than a stellar companion; earlier data showed the 
star rebrighten at mid-eclipse, altho ugh th is has weak- 
ened in recent data ijHamilton et al.ll2001|) . This indi- 
cates that the obscuring material subtends a large angle, 
which may herald the presence a low mass stellar com- 
panion to shepherd the mate rial via resonant int eraction 
or to create a distorted disk (Brvden et al. 20Q3). In ad- 
dition, the duration of the eclipse trough is lengthening 
over time, changing from 16 days to 18 days over 5 years, 
and dat a from 50-90 year s ago showed no evidence for 
eclipses l)Winn et alJ2003Tl , which means that the obscur- 
ing material is rapidly evolving. During eclipse, residual 
flux is still present at a level of ~l/20 of the flux outside 
of eclipse indicating that either (a) the obscuring mate- 
rial is porous (covering factor is < 1), (b) the eclipse is 
total, but scattered light "fills in" the eclipse, creating 
some residual flux, (c) some flux is scattered and some 
transmitted. The scattering interpretation may be fa- 
vored as the star becomes slightly bluer during eclipse as 
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one would expect from scattered light, although the fluc- 
tuations in color ar e comparable to the color difference 
ijHerbst et al.l l2003). In addition, the equivalent width of 
Ha increases dramatically from 2 A to ~ 30A indicating 
a change from an unobscured "weak-line" T Tauri star 
to a classical T Tauri star! 

The eclipse provides a natural "coronagraph" during 
which the suppression of stellar light allows any scatter- 
ing by surrounding material to become more prominent. 
We attempted to distinguish the above possibilities by 
carrying out spectropolarimetry inside and outside of the 
eclipse. Since the unpolarized flux of the star dilutes the 
percentage of polarization, we expected that the percent- 
age polarization due to scattered light would increase by 
a factor of at least 20 during the eclipse if the scatter- 
ing material lies exterior to the occulting material and if 
the starlight is intrinsically unpolarized. This scenario is 
analogous to Type 1 active galactic nuclei (AGN) which 
are unobscured and show low polarization versus Type 
2 AGN w hich are obscured and show much higher po- 
larization l|Antonucci I IT993L Some T Tauri stars show 
an anti-correlation between total flux a nd polarization 
(e.g. lMekkadeJl99ftlMenard et al.ll2003j) . which may be 
exaggerated in KH 15D as the flux variations are much 
more extreme. On the other hand, if the residual flux is 
only due to incomplete coverage by the occulting mate- 
rial, then the polarization will not change. Interstellar 
foreground polarization, which is measured to be 0.1 - 
0.4% in the V-band for NGC 2264 ijRreger Rr. Dvckll 972ft . 
will not change during eclipse. Before our observations, 
the optical polarization of KH 15D had not yet been mea- 
sured, but during eclipse it shows the characteristics of 
classical T Tauri stars which have a mean polarization of 
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Table 1. 



Table 2. 



Instrument 
& side 


A range 

(A) 


Resolution 
(A/pixel) 


Date 
CUT) 


Telescope 


Time 
(hr) 


Slit Width 


Conditions 


DBSP, red 


5620-8120 


2.45 


2003 Jan 25 


Palomar 


3 


2'.'0 


Partly cloudy 


DBSP, blue 


3350-5750 


3.38 


2003 Feb 23 


Keck 


3 


1'.'5 


Partly cloudy 


LRISp, red 


3200-5850 


2.0 


2003 Feb 24 


Keck 


1 


l'.'O 


Clear 


LPJSp, blue 


5650-9400 


2.0 


2003 Feb 25 


Keck 


2 


1"5 


Clear 



1.2% with a standar d deviation of 0.8% and maximum 
polarization of 2.5% ljMenardlll99l|) . 

In §2 we summarize the spectropolarimetric observa- 
tions, in §3 we give the results of the observations, in §4 
we discuss the implications of these results for models of 
this T Tauri star, and in §5 we review our conclusions. 

2. OBSERVATIONS 

The primary observational goal was to measure the dif- 
ference in polarization between eclipse and non-eclipse. 
Outside eclipse, the star is sufficiently bright to ob- 
serve at the Palomar 200- inch telescope (V ~ 16), but 
during the eclipse the star is a factor of ~ 20 fainter 
which required using Keck to obtain a similar number 
of photons in a reasonable observing time. The Palo- 
mar observations were ta ken with the Double Spectro- 
graph spectropolarimeter l)Oke fc Gunnlll982t iGoodrichl 
Il991|) . while the Keck I observations were taken with 
the Low Resolu t ion Imaging Sp ectrograph polarimetcr 
ijOke et al.lll995HGoodr ich 1991). Both instruments add 
dual-beam polarimetry optics to a double spectrograph 
incorporating dichroics to separate the blue and red 
beams. The observations and reductions w ere done us- 
ing st andard spectropolarimetry techniques (jMiller et al.l 
1988) . The observations of KH 15D were performed in se- 
quences of four 15-minute exposures at four different po- 
sitions of a rotating half-wave plate. Each four-exposure 
sequence allows measurement of the Stokes parameters q 
and u. For each hour of observations, the position angle 
of the spectrograph slit was aligned with the parallactic 
angle for the midpoint of the hour. The instrumental 
setup is described in Table 1, while the observational 
summary is given in Table 2. Null, polarized, and flux 
standards were observed during each night. Since instru- 
ment aj_p^larization_can be present in the Palomar 200- 
inch (Og le et al.l ll999V we observed four null standards 
each at two perpendicular slit position angles. We used 
the measured polarization of the null standards as an es- 
timate of the mean and uncertainty of the instrumental 
polarization, and we subtracted the mean instrumental 
Stokes parameters q and u from the measured Stokes 
parameters of KH 15D to correct for the instrumental 
polarization. At Palomar the instrumental polarization 
was found to be 0.30 ± 0.22% over 4600 - 5600 A and 
0.13 ± 0.13% over 5650 - 8500 A. At Keck the instru- 
mental polarization was negligible, < 0.1%. At Palomar 
the observations were obtained through clouds, while the 
Keck nights were mostly clear, but not photometric. 

A 7th magnitude B2III star, HD 47887, is located 
39" away from KH 15D. In the Keck observations a 
diffraction spike from this bright star periodically swept 
through the slit, passing within several arcseconds of KH 
15D in several exposures. In these exposures sky sub- 



Table 3. Palomar Observational Summary 



Band (A) 


p(%) 




(Fun / Fec) & 


4600-5600 


0.26±0.25 


31±27 


1.00±0.04 


6000-7000 


0.17±0.21 


-19±35 


1 


7000-8000 


0.16±0.23 


-22±41 


1.05±0.04 



a The mean of the ratio of the uneclipsed (Palomar) flux to the 
eclipsed (Keck) flux within the given band, emission lines excluded, 
divided by the ratio in the 6000-7000 band. 



Table 4. Keck Observational Summary 



Band (A) p{%) 6 P (°) 



4000-4600 
4600-5600 
6000-7000 
7000-8000 
8000-9000 



1.46±0.70 
1.63±0.19 
1.85±0.14 
2.25±0.10 
2.24±0.29 



23.9±14. 
21.8±3.3 
18.0±2.2 
16.8±1.2 
15.5±3.7 



traction was difficult, making it impossible to accurately 
measure the polarization and also adversely affecting the 
color of the extracted spectra for two of the six hours. 

Figure 1 depicts the phase of our observations with 
respect to the 48.36 day period of KH 15D. 

3. RESULTS 

The continuum polarization in the Keck and Palomar 
data was measured over wavelength bins of width 1000 A. 
We averaged the 3 measurements from Palomar and sub- 
tracted off the mean polarization measured from obser- 
vations of 4 null-standard stars. We computed the mean 
of 4 measurements from Keck and the error bars were 
computed from the standard deviation of these measure- 
ments plus the statistical error of the individual measure- 
ments. The Palomar data show a small enough polariza- 
tion and large enough errors that the data are consistent 
with being unpolarized at the ~ 1 a level (see Figure 2). 
However, the Keck data showed a much higher polar- 
ization, and are inconsistent with the Palomar data at a 
high significance. Tables 3 and 4 summarize the polariza- 
tion measurements of the eclipsed and uneclipsed data. 
The eclipsed data show a slight increase in polarization 
from the blue (4600-5600 A) to the red (7000-8000 A) 
at 3 (7 significance. Most of this change is due to an 
increase in the q Stokes parameter (Figure 2), while u 
remains constant; however, the errors are large enough 
that a change in the polarization angle is only significant 
at 1.5 a. 
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Fig. 1. — Phased lightcurve of KH 15D from 2001-2002 campaign (Herbst et al. 2002). The arrows indicate the phases of our observations. 
The solid line is a model fit to the data described in §4.3. 
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Fig. 2. — Normalized Stokes parameters q and u of KH15D. The 
squares arc the data in eclipse, while the plus signs are data out 
of eclipse. Horizontal bars represent the wavelength ranges over 
which the Stokes parameters were measured. 

Figure 3 shows the total flux spectrum in and out of 
eclipse. The spectral shape appears nearly identical in 
and out of eclipse, despite the change in total flux of 
~ 3.5 magnitudes. We measured the mean of the ratio 
of the uneclipsed to the eclipsed data in each of the three 
bands (4600-5600, 6000-7000, and 7000-8000 A), exclud- 
ing the emission lines, and computed the standard devi- 
ation of the ratio within each band, which is reported 
in Table 3. The flux ratios differ between the bluest 
(4600-5600 A) and reddest (7000-8000 A) bands by only 
1 a, which is statistically consistent with no color change. 
The absolute value of the uneclipsed to eclipsed flux ra- 
tios are not to be trusted as the nights were not photo- 
metric, so we report only the ratio of each band to the 
6000-7 000 A band. F or com parison with the photometric 
data of lHerbst etHI l)2003ft we also computed synthetic 
photometry of our spectra which gave V — R = 0.81 mag 



and V — I = 1.45 mag during eclipse and V — R = 0.79 
mag outside of eclipse, whi ch are in agre e ment with the 
range of values reported bv lHerbst et al.1 l)2003|) . We es- 
timate the uncertainty on these colors to be ~ 0.1 mag as 
there appear to be differences between individual spectra 
on scales of ~ 500 A which may be due to fluctuations 
in the flat-fields, difficulty of sky subtraction, overall sys- 
tematic uncertainties in flux calibration, or other instru- 
mental variations. 

There are a large number of stellar photospheric fea- 
tures which are seen in both spectra and are indicative 
of a K6V spectral type with T e ff ~ 43 00 K, somewhat 
hotter than the K7V type estimated bv lHamilton et all 
l)2001[) which may be due to the fact that our spectrum 
covers bluer wavelengths. Almost all of these features 
show the same equivalent widths in both spectra, indicat- 
ing that the eclipsed spectrum is a scaled version of the 
uneclipsed spectrum. We hav e measured the Lick indices 
l)Worthev fc Ottavianil 119971) for both the eclipsed and 
uneclipsed spectra, plotted in Figure 4. The H<5 indices 
are contaminated by the [SII] emission line in the Keck 
data. The Na D equivalent width is smaller during the 
eclipse, which we attribute to unobscured Na D emission 
filling in the absorption line. The measured indices for 
KH 15D agree best with those of a K6V star; the agree- 
ment is good enough that there is no evidence for optical 
veiling. The Ha emission line has a much higher flux in 
the eclipsed spectrum than in the uneclipsed spectrum by 
a factor of ~ 3.3 (since our data were not photometric, 
this factor is probably ~ 2.7). This may be partly due 
to vari ability and partly due to partial eclipse of the Ha 
line as lHamilton et al.l l)2003|) found that the flux of Ha 
could decrease or increase during eclipse and the shape 
of Ha varied significantly as well. They measured an 
equivalent width of 2 A outside eclipse and 40-70 A in 
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Table 5. Emission Lines 



Line 


Flux a 


EW 


(A) 


(10 -15 erg cm -2 s _1 ) 


(A) 


Ha A 6563 b 


5.88 ±0.71 


3.2±0.39 


Ha A 6563 


1.74 ±0.19 


29.7±3.2 


[S II] A 4069 


0.362±0.025 


22.2±1.5 


Na D A 5890 


0.153±0.015 


3.6±0.36 


[S II] AA 6716,6731 


0.093±0.013 


1.6±0.23 


[O I] A 6300 


0.397±0.065 


7.0±1.1 


[0 I] A 6364 


0.132±0.017 


2.3±0.3 



a The absolute value of these measured fluxes is unreliable due to 
non-photometric conditions and slit losses. 
b Uneclipsed line (the rest are eclipsed). 



eclipse, compared with our measurements of 3 A and 30 
A, respectively. In addition, four emission lines appear 
in the eclipsed spectrum that are not apparent in the 
uneclipsed spectrum. The fluxes and equivalent widths 
of the emission lines are given in Table 5, which we have 
computed by subtracting off a continuum template, with 
errors estimated from the standard deviation of the resid- 
uals. Included is the Na D absorption line at A5890 A 
which is stronger in the uneclipsed spectrum than in the 
eclipsed spectrum (we have used the uneclipsed spectrum 
as the continuum template for this line). 

4. DISCUSSION 

The variation of the polarization between eclipse and 
non-eclipse indicates that the scattering region is not 
completely obscured by the occulting material. If the oc- 
culter completely blocks all the light from the star, if the 
scattering is nearly isotropic, and if the scattering region 
is nearly spherically symmetric, then the residual flux, 
~ 5%, indicates that the optical depth of the unobscured 
scattering material is ~ 0.05. Since the spectra inside 
and outside eclipse are consistent with no change in the 
spectral shape, and since the spectrum within eclipse is 
consistent with very weak wavelength dependence of the 
polarization, we conclude that the scattering mechanism 
is nearly achromatic, implying that the scatterers may 
be electrons or large grains. We discuss each of these 
possibilities in turn. 

4.1. Electron scattering 

If the scattering mechanism were electrons, then the 
low measured polarization of 2% indicates that the scat- 
tering region need be only slightly asymmetric . For ex- 
ample, in the model of lBrown fc McLeanl 1)19771) . the po- 
larization caused by an axisymmetric Thomson scatter- 
ing envelope is P = sin 2 i/(2a + sin 2 i), where a is a 
"shape factor" which depends on the precise geometry 
of the scatterers. Since the eclipsing region is likely in 
the equatorial plane of the circumstellar material, we can 
assume that i ~ 90°. Thus, a ~ 25 for P ~ 2%. For 
a scattering geometry which consists of a slightly oblate 
or prolate spheroid of axis ratio a = 1 + e, then a ~ 5/e, 
so |e| ~ 0.2 to create the observed polarization, i.e. only 
a 20% deviation from spherical symmetry is required to 
give the measured polarization. Note that this example 
is simply illustrative since we have no preferred model 
for the scattering region. 



If we assume that all of the light during eclipse is due to 
electron scattering, then we can estimate the bound-free 
continuum emission expected from this gas. Consider a 
nearly symmetric region of uniform electron density and 
temperature surrounding the star of size r es . We assume 
that the eclipsing material is located near the Keplerian 
radius, rj^ep, for a 48 day period. The low polarization 
implies that the electron scattering region must be at a 
larger radius, r es > rx e p ~ 0.2 AU to avoid becoming 
too asymmetric due to obscuration. Since the material is 
nearly symmetric, the mean optical depth of the scatter- 
ing region is r es ~ 0.05, which implies an electron num- 
ber density of n e < 2 x 10 10 cm~ 3 (assuming a uniform 
density). The bound- free emission expected at, say, 6000 
A must be less than F v ~ 10~ 28 erg cm~ 2 s~ 1 Hz~ 1 to 
avoid changing the strength of the Lick indices. Comput- 
ing the emission from an optically-thin spherical region 
of plasma in LTE, we find that T > 2 x 10 6 K to be con- 
sistent with this upper limit. Note that this is a strong 
lower limit since it is likely that the scattering region 
need be larger to give the small observed polarization. In 
addition, one can place a constraint on the temperature 
based on the strength of Balmer line emission; however, 
the interpretation is complicated by the unknown opti- 
cal depth of the Balmer lines. Even so, the lower limit 
on the temperature based on the lack of detected bound- 
free emission is four times the virial temperature at rxept 
meaning that this gas would be strongly unbound and 
will outflow in less than a dynamical timescale, implying 
a mass loss rate of 1O _7 M0 yr -1 , which is rather large 
for a weak-line T Tauri star. Thus, it is more likely that 
the scattering is not due to electrons, but may be due to 
large dust grains, which we discuss next. 

4.2. Dust scattering 

If the scattering were due to dust, then the grain sizes 
would have to be fairly large to allow the weak wave- 
length dependence of the polarization and flux. This is 
true for the zodiacal light which reflects a n optical spec- 
trum that is identical to that of the Sun ijLeinert et alJ 
1998). Based on the achromaticity of the zodiacal scat- 
tering, the scatterers are thought to be large dust grains, 
10-100 /im in size, and there is further evidence that these 
grains might be quite "fluffy" (i.e. low filling factor) with 
densities as low as pf ~ 0.25 g cm -3 . 

For the estimation of the wavelength dependence 
of the scattering and polarization by dust particles 
we approximate the particles as spheres (Mie scat- 
tering) and compute the polarized differential scat- 
tering cross sec tion by single partic l es using a stan- 
dard approach iVoshch innikov) 120021 iWiscombel 119801 
iBohren fc HufFmanl 11983 iWolflboOpfT 6 We have chosen 
the distribution of refractive i ndices of the materials from 
iWeingartner fc Draind 1)2001)) : 62.5% astronomical sili- 
cate and 37.5% graphite. Since graphite is crystalline, it 
was considered as two "different" materials in the com- 
putation: the first with a refractive index as measured 
perpendicular to the crystal rotation axis, 25%, and the 
second parallel, 12.5%. Wc assumed that distribution of 

6 The code for calculation of the Mie scattering coefficients is 

available at 

http://mc.caltech.edu/~swolf/miex-web/miex.htni - or 
contact the authors . 
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Fig. 3. — Spectrum out of eclipse (top line) and in eclipse (bottom line). 
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Index 

Fig. 4. — Lick spectral indices for the eclipsed (blue crosses) and uneclipsed (red squares) data, as well as a K7V star (solid line), K6V 
star (dotted line), and K5V star (dashed line). Note that the molecular band indi ces with units of magnitu des (CN, Mg, and TiO) have 
been multiplied by a factor of 10 for clarity. The index numbers are taken from Worthcv & Ottaviani (1997). 



the dust size, a, follows a power-law dn/da oc a a for ties o f the ensemble of particle species and sizes ijMartirJ 
a m in < a < CLmax, and derived the scattering proper- 119781 Wolf & Voshchinnikov, in prep.). 
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We averaged the scattered polarized flux over the entire 
47r steradians divided by the total scattered flux to esti- 
mate the mean polarization as a function of wavelength 
for each dust model. We then compared this to the rela- 
tive wavelength dependence of the measured polarization 
(the absolute polarization depends on the exact geometry 
of the scattering region, which is unknown; see the next 
section for one particular model) . We also compared the 
relative wavelength dependence of the scattered flux by 
dividing the eclipsed spectrum by the uneclipsed spec- 
trum. Given the assumptions that (a) the wavelength 
dependence of the flux and polarization does not depend 
on the geometry of the scatterers, and (b) the scatter- 
ers can be described by a single power-law, we can then 
place a constraint on the parameters of the dust model. 
We find that the lower limit, a min , does not change the 
fit much, so we simply vary l/im< a max < 100/im and 
< a < 2. The minimum x 2 is 12 for 5 data points 
for a max = 8/Ltm and a = 0, with the confidence con- 
tours plotted in Figure 5. At 90% confidence, a < 2 
and bfim< a max < 10/im, which indicates that the scat- 
tering is dominated by large dust grains. Unfortunately, 
the minimum reduced Xr — 2-2, which is rather large. 
This is due to the fact that the measured polarization 
increases toward longer wavelengths, while the dust scat- 
tering models predict a polarization that is nearly con- 
stant, or slightly decreasing, with wavelength. 



100 



i 10 




2.0 



Fig. 5. — Confidence limits for dust models as a function of a ma x 
and a. 



The opacity of large dust grains scales as n/p = 
0.75a _1 p _1 ~ 3 x 10 3 a^ 1 (pf / p) cm 2 g -1 , where a w = 
a/10/im. Thus, the Eddington limit for the smaller of the 
large dust grains scales as ^Edd, dust = l-7xl0 34 aio(p//o/) 
erg s _1 , which is similar to the luminosity of KH 15D, 
£khi5D ~ 2.4 x 10 33 erg s _1 . Thus, smaller dust grains 
will be removed from the system by radiation forces. In 
addition, radiation forces on the larger grains may be 
significant enough to modify the structure of the system. 

4.3. Scattering Model 

Hcr bst et al.l l)2003l) have proposed a model for the oc- 
culter which consists of a semi-infinite edge crossing the 
star. This model is attractive as the finite size of the 
star and limb darkening can account for the shape of the 
ingress and egress, and, if the occulter is at the Kcple- 
rian radius, then it must be inclined to the orbital plane 
by 15° to give the correct duration of ingress and egress. 



However, we have found that the ingress and egress of 
the eclipse (Figure 1, near phase of 0.3 and 0.7) show 
a more gradual rolloff than expected for a semi-infinite 
occulter. This can be fit with an optical depth of the oc- 
culter which scales with position as r cx x~ 2 , where x is 
the coordinate perpendicular to the edge of the occulter, 
which gives an excellent fit to the phased lightcurve data 
(solid line in Figure 1). 

This material causing extinction will scatter as well, 
and in the case of large grains there will be significant 
forward scattering (van dc Hulst 1981). In fact, this may 
account for the scattered flux present during the eclipse; 
however, to compute the scattered light requires a knowl- 
edge of the three-dimensional distribution of scattering 
and occulting material, as well as a knowledge of the 
scattering properties. We construct such a model by as- 
suming that (a) the occulter is located at the Keplerian 
radius and (b) the 15° tilt of the occulter is due to a 
warped disk. We model the material as a succession of 
tilted circular annuli, creating a warped disk, with the 
tilt angle given by 

(r - r w ) 2 ~ 



Hilt 



, exp 



2al 



(1) 



where r w is the radius of maximum tilt, i w is the maxi- 
mum tilt, r is the radius in the disk, and o~ w is the width 
of the tilted region, and we assume that the nodes of 
the warped annuli are aligned. The observer we place at 
an angle i b s to the equatorial plane, and we rotate the 
warped disk as a solid body with a period equal to that of 
KH 15D. In addition, we cut the disk off at r,„ = 0.5r w 
and r ou t — 4r tlI . The opacity of the scattering mate- 
rial we distribute as k — Ko[h(r /r w ) / (6 — n + itut)] 2 , to 
produce the observed roll-off in intensity, where h is the 
angular scale height, 9 is the polar angle, and kq is a 
fiducial opacity. We find a good fit to the lightcurve for 
i w = 17.6°, i obs = 6.7°, h = 1.6°, r w = 0.22 AU (this 
is the Keplerian radius), a w — r w . When taking into 
account only extinction, this distribution produces the 
solid lightcurve in Figure 1 after the addition of 4.5% of 
the uneclipsed flux during eclipse. However, part of the 
extinction is due to scattering, so we also have computed 
the singly-scattered component of this density distribu- 
tion assuming an albedo of unity and using the sca ttering 
function for zodiacal light derived by Hong (1985) which 
is strongly forward scattering and has a maximum po- 
larization of 22%. We find that the polarized flux of 
the scattered light during eclipse is 0.1% of the unob- 
scured total flux, which is the same as the observed po- 
larized flux during eclipse of KH 15D (Figure 6). Figure 7 
shows the appearance of the scattered light in this model 
at mid-eclipse (an animated version of this figure which 
shows the appearance as a function of phase is available 
on-line) . 

One prediction of this model is that the polarization 
angle should vary by ~ 15° during the eclipse. How- 
ever, the total model flux during mid-eclipse is only 2.5% 
of the uneclipsed flux which is about half the observed 
mid-eclipse flux, indicating that there may be another 
component of scattered light which has low polarization, 
such as diffrac ted light, wh ich is not taken into account 
in the model of Hong (1985), or back-scattered light. In a 
future study we may carry out more sophisticated mod- 
eling, varying the radius of the warp as well as the sym- 
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Fig. 6. — (a) Flux of KH 15D as a function of phase from lHerbst et a l. 1 2003), along with our model (solid curve) and the singly-scattered 
light (dashed curve), (b) Model polarized flux normalized by Fo at 7000-8000 A (unity would be 100% polarized light). A t mid-eclipse, 
we ha ve taken our observed polarization and multiplied by the mean flux of the I-band observations near mid-eclipse by IHerbst et alj 
(2003), F/Fq = 0.047 ± 0.01, since our observations were not photometric. The crosses with error bars show our observations, which 
agree remarkably well with the model, (c) Model polarization. Note the disagreement at mid-eclipse which may be due to an unpolarized 
component which also produces the disagreement in the total flux at mid-eclipse, (d) Polarization angle. We have not plotted the observed 
angles since we do not know the plane of the disk on the sky and our temporal coverage is not sufficient to determine the shape of 0(t). 
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metry of the warp, and including multiple scattering and 
varying the dust grain model, to see how the scattered 
flux varies with dust composition and with wavelength 
of the observed light, and we will see whether diffrac- 
tion, multiple-scattering, or back-scattering can explain 
the brightening at mid eclipse. There is evidence that 
the ingress and egress differ slightly in time between the 
even and odd eclipses, which may imply that the warp is 
an m = 2 mode rather than to = 1 as in our model. If so, 
the back warp will be partly eclipsed by the front warp, 
while the edges will be visible during mid-eclipse, which 
may explain the rebrightening. Alternatively, the sym- 
metry during mid-eclipse of the front warp may cause a 
brightening due to diffraction, which is likely to be wave- 
length dependent. If the albedo is less than unity, then 
mid-infrared observations may point toward an m = 2 
configuration if they also show a peak at mid-eclipse since 
the infrared light is emitted isotropically, similar to the 
back-scattered light. 




Fig. 7. — Scattered light at mid-eclipse for the model described 
in the text. 



5. CONCLUSIONS 
A summary of the results: 

f) Near mid-eclipse the polarization is ~2%, and is 
consistent with zero outside of eclipse (0.2±0.2%). Com- 
bined with other photometric data, the polarized flux 
near mid-eclipse is ~ 0.1% of the flux outside of eclipse. 
This indicates that there is an unobscured scattering re- 
gion surrounding the T Tauri star which creates the po- 
larization; the polarization cannot be solely due to fore- 
ground transmission effects (not associated with the star) 
as there would no variation. 

2) The data are nearly consistent with no wavelength 
dependence of the polarization during eclipse. This in- 
dicates that the scatterers are approximately grey. In 
addition, electron scattering is unlikely to cause the ob- 
served polarization. 

3) The spectrum within the eclipse is consistent with 
being a scaled version of the spectrum outside of eclipse 
with no wavelength dependence of the ratio except for the 



presence of stronger emission lines which also argues for 
grey scattering; Rayleigh scattering is ruled out. Our Mie 
dust scattering model implies a < 2 and 5/im< a max < 
10/im at 90% confidence; however, the best fit model is 
not an adequate fit. 

4) The qualitative features of the eclipse can be re- 
produced by a model of a warped disk that orbits in 
front of the star, periodically eclipsing the star, having 
an extended atmosphere which causes a decrease in light 
before and after eclipse as well as scattered light within 
eclipse, which is consistent with the observed polarized 
flux. 

The polarization adds an extra clue that the eclipsing 
material completely occults the star since the maximum 
observed polar ization of classical T Tauri stars is ~ 2.5% 
l|Menardl 19911) and the KH 15D spectrum looks like that 
of a classical T Tauri star during eclipse. This adds to 
the evidence that the eclipse is complete based on (a) the 
steep egress and ingress; (b) the flat bottom of the trough 
of the eclipse; (c) the fine tuning that would be required 
to arrange 95% covering. It is interesting that the trough 
of the eclipse does show some slight variations, namely 
an increase near mid-eclipse, which may indicate that 
the scattering region is partially eclipsed and this eclipse 
varies with time. It would be interesting to measure the 
polarization as a function of time within the eclipse to 
look for variation. 

We have modeled the eclipse with a warped disk since 
a discrete cloud would be sheared in ten orbits, while the 
eclipses have repeated with similar durations for ~ 50 or- 
bits. However, this does not mean that the warped disk 
model is without its own problems: If the warped model 
is correct, what drives the warp at such a short period? 
The interaction of the stellar magnetosphere is proposed 
to cause a similar occulting beh avior in the T Tauri star 
AA Tauri (Bouvic r et al.lll999lk however, this star has 
a much short period between occultations, 8 days, while 
the longer period of KH 15D is much longer than the spin 
period of any known T Tauri star. One can place a lower 
limit on the warp radius of r m i n = 10i?* since inwards of 
this radius the warp would have to be at > 45° to create 
the correct duration of ingress due to the finite size of the 
star. At r m i n , the Keplerian period is 0.8 days, while a 
warp precession timescale is typically much longer than 
the Keplerian times cale. It may be tha t a planet excites 
a wave m the disk <|Winn et all 12003(1 o r that an elon- 
gated dusty vortex causes obscuration l|Barge fc Vitonl 
2003); each of these models might be consistent with the 
observed polarization profile, which remains to be quan- 
tified with a scattering calculation similar to the one pre- 
sented here. 

Future observations of interest would be to (1) mea- 
sure the magnitude and angle of polarization throughout 
ingress and/or egress to test the warp model; (2) mea- 
sure the polarization at higher spectral resolution within 
the Ha emission line to see whether it is unpolarizcd; 
(3) look for variations of the flux during eclipse to de- 
termine whether the stellar variability is washed out by 
the time delays from scattering (although this requires 
a rather extended scattering region). In addition, near- 
and mid-infrared observations may show fluctuations if 
the occulting region has a warped geometry, as we have 
suggested. 
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